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Abstract: Thermodynamic properties of non-relativistic bosons and fermions in two spatial 
dimensions and without interactions are derived. All the virial coefficients are the same except for 
the second, for which the signs are opposite. This results in the same specific heat for the two 
gases. Existing equations of state for the free anyon gas are also discussed and shown to break 
down at low temperatures or high densities. 



1 Introduction 

Physics of two-dimensional systems used to be of primarily academic interest in 
providing mathematically simpler versions of realistic problems in three dimensions. 
The fundamental discovery of Leinaas and Myrheim in 1977 of the possibility for 
intermediate quantum statistics for identical particles in two dimensions interpolat- 
ing between standard Bose-Einstein and Fermi-Dirac statistics therefore generated 
initially little general interest. Particles with this new statistics were later named 
anyons by Wilczek 0. He showed that they could be modelled as carrying both a 
charge and a magnetic flux locked together in a very special way. 

With the progress of modern microelectronics this situation changed 0. The 
discovery of the integer quantum Hall effect @] could be explained by the special 
quantum effects of electrons effectively confined between layers of different mate- 
rials in planar transistors. Later it became clear that excitations in the fractional 
quantum Hall effect obeyed intermediate statistics 0. This in turn suggested that 
the mechanism behind the new, high-temperature superconductivity seen in layers 
of copperoxide planes in different materials, could also be anyonic ]§, 0, 0- 
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Even if anyons should not turn out to be of great practical importance in the 
coming years, two-dimensional systems in general will play an important role in 
microelectronics and material sciences. Since the third dimension is effectively frozen 
out only at very low temperatures, it is also obvious that quantum effects will govern 
the behaviour of these new systems. 

Here we consider the thermodynamics of ideal quantum gases in two dimensions. 
Ordinary bosons and fermions are treated in the next section where their equation 
of state is derived. In Section 3 we review the basic properties of anyons and solve 
the simplest case of two anyons bound by a harmonic potential. The resulting 
energy spectrum is then used in Section 4 to calculate the second virial coefficient 
for anyons. Including some recent results for the higher virial coefficients, we then 
discuss in the last section their full equation of state. 



2 Statistical mechanics of bosons and fermions 



We consider N free and identical particles of mass m in a box of two-dimensional 
volume V. Each particle is characterized by the momentum k and corresponding 
energy e = k^/2m. When they are in thermal equilibrium at temperature T and 
chemical potential n, their thermodynamics follows from the grand canonical parti- 
tion function 

E = l[[l±e-^^''^-^'^f\ (1) 

k 

where the upper signs are for Fermi-Dirac statistics and the lower ones for Bose- 
Einstein statistics 0. Here (3 = l/ksT where fc^ is the Boltzmann constant. In this 
ensemble the corresponding free energy density is simply the pressure P = InS/ j3V 
which becomes 

/3P = ±:^^ln[l±e-^(")]. (2) 

k 

The particle density p = N/V then follows from 
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Similarly, if E is the total energy of the particles, the energy density S = E/V is 
given by the derivative 

Since the lowest one-particle energy is here e = 0, the chemical potential can not 
become negative for bosons. Their fugacity z = exp will thus always be less 
than one. No such upper limit exits for fermions. 



In two dimensions we can replace the momentum sum by the integral 



in the thermodynamic limit where V oo. Writing kdk = mde, we see that the 
particle density (j^) is given by the integral 



2nh^ Jo e/5(^-/^)±l 27r/3r 

Introducing the thermal wavelength A = (27r/5/;,^/m)^/^, we can write the result for 
the fugacity as 



^ = ^ ^1 _ e^^^^j . (7) 

The upper signs are still for fermions and the lower ones for bosons. For bosons 
z smoothly approaches the value one as the temperature goes to zero. Thus there 
is no Bose-Einstein condensation in two dimensions except at T = when all the 
particles are in the ground state. 

For the energy density (^) we similarly find in the thermodynamic limit 

S = ^ de--rA • (8) 

A2 Jo e/3(^-/^) ±1 ^ ^ 

A partial integration where the boundary term vanishes then gives 

1 roo 

£ = ±— de\n[l±e^^^'~''^]. (9) 

This is seen to equal the pressure (|^) in the same limit. We thus have P = S. In three 
dimensions the corresponding result is the more well-known relation P = {2/3)S for 
non-relativistic particles. 
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The integral for the pressure can be expressed in a more compact way. 
Changing variable of integration to t = 1 ± 2;exp we obtain 

/3P(A/i) = T^Li2(l±^) (10) 



where Li2(a;) is the dilogarithmic function defined by ||T0 



log t 

U,{x) = - dt-^. (11) 
Ji t — I 

When < X < 2 one can obtain a convergent series for the pressure by expanding 
the logarithm in (|^) and integrating term by term. This expansion is thus valid at 
all temperatures for bosons and only at high temperatures for fermions. 

With the explicit solution (^ for the fugacity we can now obtain a differential 
equation for the pressure which will give the equation of state. Forming the partial 
derivative 

we see that the first derivative on the right-hand side is just the density (0) and the 
second is obtained from (^ as 



The equation of state is thus explicitly given by the integral 



PP = tL dp'-^^^, (12) 



as first obtained by Sen |TT|. 



At high temperatures we can express the equation of state in terms of the virial 
coefficients. They can now be obtained by using the expansion 

-^ = T.Bn^ (13) 



which is convergent for < |x| < 2n ||T0|. It defines the Bernoulli numbers Bn, 
giving Bo = 1, 5i = -1/2, B2 = 1/6, Bs = 0, ^4 = -1/30 etc. with B2n+i = 0. 
We then obtain from ( [T^ ) the virial expansion 

00 

I3P = Y.A,/A'^'-'^ (14) 



4 



where the dimensionsless virial coefficients are simply given in terms of the Bernoulh 
numbers |Tl|, |T^ , 



A, = (tY-'^. (15) 

The lowest ones are A2 = ±1/4, A3 = 1/36, = 0, etc. and give for the equation 
of state at high temperatures or low densities 

/3P = p± -p'A' + -p=^A^ - -^p'A^ + ^^p7Ai2 + . . . . (16) 
^ V ^ Se'^ 3600'^ 211680'^ ^ ' 

Notice the interesting fact that since all the odd Bernoulli numbers are zero after Bi, 
all the even virial coefficients are zero after A2. Hence, the only difference between 
the pressures in the bosonic or fermionic gases comes from the second term in the 
expansion (|16|), i.e. 

Pf-Pb = ^p'A' = -V. (17) 

The result does not depend on the convergence of the virial expansion. It follows 
directly from (|T0|) when we make use of the special property of the dilogarithmic 
function|ll0[ 



Li2(x) + Li2(x"^) = -llog^x, (18) 

where x = exp (pA^) in our case. The difference is just the pressure or energy density 
£0 = fi^Tip^ /m in the Fermi gas at zero temperature as shown in the Appendix. Since 
it is constant with respect to T, we see that the specific heats of these ideal quantum 
gases are equal as pointed out by Aldrovandi |TB[. It would be of interest to have a 
more direct understanding of this simple result. 

We have plotted the pressures of the gases as a function of temperature at fixed 
density in Fig.l. At high temperatures we can use the virial expansion (p!^). Since 
it only converges for pA^ < 27r, it breaks down at low temperatures. For the Fermi 
gas just above T = we can instead use the Sommerfeld expansion ||^ derived in 
the Appendix. In two dimensions it is found to terminate after the second order. At 
slightly higher temperatures it smoothly matches up with the virial expansion. We 
can obtain the pressure of the Bose gas by just subtracting the constant term £q from 
the Fermi pressure. In the Appendix we have also derived exact equations of state 
for these g explicit functions of the dimensionless quantity pA^ combining the 

results in (0) and (|^. 
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3 Quantum mechanics of anyons 



The wavefunctions of bosons and fermions differ under the interchange of two or 
more particles. While the wavefunction for bosons is completely symmetric, the 
wavefunction for fermions is completely antisymmetric. An interchange of two par- 
ticles described by the wavefunction '?/'(ri,r2) is shown in Fig. 2a. Their relative 
vector is rotated by the angle A0 = vr. The wavefunction is then transformed into 
'?/'(r2,ri). Since the two particles are indistinguishable, this wavefunction can only 
differ by a phase angle 9 from the initial wavefunction, 

V'(r2,ri) = e'V(ri,r2). (19) 

Had the interchange taken place in the opposite direction as in Fig. 2b, the final 
wavefunction would have been 

V^'(r2,ri) = e-^V(ri,r2) (20) 

since these two exchanges followed by each other are equivalent to no exchange. 

In three dimensions we now see that this latter interchange is equivalent to the 
first since it can be obtained by an out-of-page rotation around the line connecting 
the two particles in their final position. Thus we have '?/''(r2,ri) = '?/'(r2,ri) which 
implies that exp(2i^) = 1. This equation has the two solutions = mod 2tt 
and 6 = 7T mod 2ti. In the first case the wavefunction is symmetric under the 
interchange and the particles are bosons, while in the other case it is antisymmetric 
and the particles are fermions. 

In two dimensions the two interchanges cannot be related for particles and the 
parameter 9 can have any value in the interval < < 27r. We then have the 
possibility for intermediate statistics which interpolates between Bose-Einstein and 
Fermi-Dirac statistics. The wavefunction for the corresponding particles which are 
called anyons, will thus be neither completely symmetric nor antisymmetric. Need- 
less to say, the statistical mechanics of anyon gases will be much different from the 



behaviour of bosons and fermions 14, 15 



A quantum mechanical description of anyons can be obtained by modifying 
the classical Lagrangian 



1 ^ 



L = -5:mr2-r(ri,r2,---,r;v) (21) 



1=1 



where the last term describes the potential energy of the particles and is symmetric 
in all its arguments. The quantum mechanical motion of this system is now given 
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by the integration of exp {i J dtL/h) over all possible paths the system can take pO . 
When the motion corresponds just to an interchange of two particles, the net effect 
should be a change of the wavefunction by a phase as in (0). This effect can now 
be built into the above Lagrangian by adding a topological term 

Le = L + -hY,^ij. (22) 

This additional term is the sum of time derivatives of the relative angles of all pairs 
of particles in the system as illustrated in Fig. 3. These angles obviously depend 
on the positions of the particles. When the angle 4>ij changes by Acpij = n, the 
corresponding change in the wavefunction is 

as desired. From this modified Lagrangian we can obtain the canonical momenta 
and quantization can be performed by standard methods 



In this way we can describe anyons as ordinary particles in two dimensions 
with a special, topological interaction. Bosons obey symmetric statistics and prefer 
to be in the same quantum state. Fermions can be considered as bosons with a 
statistical interaction of strength a = d/ir = 1. This results in an effective repulsion 
between fermions which is usually explained by the Pauli principle. Similarly, one 
can describe bosons as fermions with an additional, statistical interaction with a = 
1. The effect of this topological interaction for arbitrary values of the parameter a 
in systems of many anyons is not yet known. 

Essentially the only known system where exact results can be derived is the case 
of two anyons in a harmonic oscillator potential 

\/(ri,r2) = |mo;2(r? + r2). (23) 

This is equivalent to the case of two anyons in a uniform magnetic field, see e.g. [|T6|. 
Introducing the center-of-mass coordinate R = (ri + r2) and the relative coordinate 
r = ri — r2 = (r cos 0, r sin 0), the Lagrangian (P^) becomes 

Le = |MR2 + iMu;2 ^ i^^-2 ^ ^2^2 ^ ^2^2^ ^ ^^^^ ^24) 

We have here introduced the total mass M = 2m and the reduced mass fi = m/2. 
The statistical parameter enters only in the relative motion described by the last 
term in the Lagrangian. The center-of-mass motion is described by the first part 
which is just an ordinary, two-dimensional harmonic oscillator. In order to quantize 
the system we need the canonical momentum pr = mr and 

dL 

Pff) = — - = fxr'^cj) + ah. (25) 
(90 
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They allow the construction of the Hamiltonian function of the system. In the 
quantum description the relative motion is governed by the operator 



H = -- 



2/x 



I d 1 



la 



(26) 



The Schrodinger equation is separable, and the corresponding wave functions can 
be written as 



V'(r,0) 



R{r) 



(27) 



where the angular quantum number m is an even number if the particles with the 
statistical interaction are bosons and an odd number if they are fermions. For a 
given angular momentum the radial Schrodinger equation then simplifies to 




(f Id 1 , .n 



+ l/icu^r^j R{r) = ER{r) 



{2t 



which has exactly the same form as for an ordinary, two-dimensional harmonic 
oscillator with a angular momentum i = m — a. The eigenvalues are thus PI 



Enm = huj[2n + \m — a\ + 1], 



(29) 



where n = 0, 1, 2, ■ ■ ■ is the radial quantum number. If the particles are bosons, the 
angular momentum takes the values m = 0, ±2, ±4, • • •. The quantized energy levels 
as functions of the statistical parameter are shown in Fig.4. They are seen to repeat 
themselves for other values of a with period 2 although the energy of each indiviual 
state with definite quantum numbers (n, m) is not periodic. 

From Fig.4 we see that the energy levels fall into two separate classes with 
opposite slopes. Within each class the degeneracy increases linearly with the number 
of the level counted from below. For more than two anyons we only know the lowest 
energy levels from numerical calculations [|17l. Again one finds that they can be 



grouped into different classes depending on their slopes as functions of a. This can 
be explained by crude, semi-classical considerations |T^. In the case of anyons, 
exact expressions for the ground state and some excited states have been found, see 
for example [|T6|. Apart from that, few exact results are known for systems of more 
than two anyons. 
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4 The second virial coefficient 



Gases of fermions or bosons at very liigli temperatures beliave as classical, ideal 
gases, independent of the quantum statistics. We expect to see the same behaviour 
in a gas of anyons. At lower temperatures quantum effects come into play depending 
on the statistics. They result in a non-zero value for the second virial coefficient, 
which for a gas confined to a volume V, can be calculated from |^ 



B2= lim -1-2-1 (30) 




where Zi and and Z2 are the partition functions of one and two particles respectively. 
While Zi is independent of statistics and thus the same for all free particles, Z2 for 
anyons is more difficult to derive. It was first calculated by Arovas et al. [|l^ using 



Feynman's path integral formalism ||20|. Since we know all the energy levels ( pQ]) of 



the two-anyon system, it can also be directly obtained from 



nm n=0 s=-oo 

where we have written the angular momentum m = 2s when the anyons are de- 
scribed in terms of bosons. The factor Zi in front gives the contribution from the 
center-of-mass motion. From (p^) we see that it equals the partition function for 
one particle in a two-dimensional harmonic oscillator potential, i.e. 

/ 00 1 \ ^ 1 
= y e-'^^"("+2) = 27 TT- (32) 

In ( PH ) the double sum is just the product of two such geometric series and gives 
similarly 

cosh/5/i^(l - a) 
- 2sinh2(/5;.c.) ^^^^ 

for two anyons in a harmonic oscillator potential. 



The infinite volume limit ^ 00 in ( pO| ) is replaced by the zero slope limit 
a; — > for particles confined instead by a harmonic potential. In this limit (|32D 
approaches Z\ l/{f3hujY which should be compared with the result for a free 
particle in a box of volume V , 
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where A again is the thermal wavelength. Comparing these two results, we have the 
relation 



V 



' A ^ 



(35) 



In addition, in the standard formula PD| ) we must make the replacement V/2 ^ V 
when the particles are confined by an oscillator potential as here p2|. Expanding 
the last term to second order in uj, we obtain 



2^ 



The second virial coefficient for anyons can thus be written as B2 = A2A'^ where its 
dimensionless value is 



A, 



1 _ 1(1 

4 2 v-^ 



a 



(36) 



For ordinary bosons which have a = 0, we recover the result A2 = —1/4 from 
Section 2, while a = 1 gives the standard fermion result A2 = 1/4. A2 is plotted 
as a function of the statistical parameter in Fig. 5a where it is seen to interpolate 
smoothly between the fermionic and bosonic values. The cusps at the bosonic points 
are found to disappear when there is an additional interaction potential between the 
anyons depending on the separation as 1/r^ [El]. 



5 Equation of state for anyons 

In order to calculate the higher virial coefficients B^, B4, . . ., we need the partition 
functions for = 3,4,... anyons. But the energy levels and the corresponding 
degeneracies for these systems are not known. Instead of such exact calculations 
one has to rely on perturbation theory or numerical methods. Following the latter 
approach, Myrheim and Olaussen |2^ have used Monte Carlo simulations of the 
three-anyon problem which have enabled them to extract the third virial coefficient. 
In this way they have obtained the simple result 
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Corrections are expected to go as sin^ 6, but are absent or so small that they cannot 
be determined with the present level of accuracy in the simulations. Higher virial 
coefficients will be very difficult to obtain by this method. ^3 is plotted in Fig. 5b 
and is seen to take the common value = 1/36 derived previously for fermions 
and bosons. 

The Monte-Carlo result is in agreement with analytic calculations where the 



grand partition function has been obtained perturbatively to ffist |25] and second |26 



order in the statistical angle. These results have more recently been conffimed in the 



Chern-Simons formulation of anyons ||T5[ to the same order in perturbation theory 



2% . Extracting the lowest virial coefficients and plotting the resulting pressure as 



function of temperature at fixed density, one finds as expected a curve in between 
the fermion and boson lines in Fig.l as long as pA^ < 1 [0]. It has been shown 
that the relation S = P between energy density and pressure is valid for anyons as 
it is for bosons and fermions [^. One can thus also obtain the specific heat for 
anyons at high temperatures [^]. Since the virial expansion cannot be used at low 
temperatures, one obtains no information about the thermodynamics of anyons in 
this regime where quantum effects dominate. 

Instead of relying on the virial expansion, one can try to use the pressure as 
function of fugacity in the discussion of the equation of state. To ffist order in 9, 
this is found to be p5l E7 



Lioil 



TX 



{31 



With the density p = zd/dz{pP) we can plot the pressure as function of tempera- 
ture with the fugacity as a parameter. The result is shown in Fig. 6 for two small 
values of the statistical angle. We see that when the temperature gets sufficiently 
small, the resulting pressure is no longer well-defined [^]. When 6 becomes smaller, 
this happens at a corresponding lower temperatures. A similar breakdown of pertur- 
bation theory is also seen using the analytical results |]26| , ^ for the pressure valid 
to second order in 6. It can be traced back to a violation of dp/dz > which follows 
from general statistical mechanics for systems in thermodynamic equilibrium p5 



This inequality is derived in the Appendix. Since it is here perturbatively violated, 
one can even get negative particle densities for sufficiently low temperatures. Higher 
order perturbative results will probably not give much more understanding of the 
physical phenomena causing this unexpected behaviour at low temperatures in the 
anyon gas. 

When the statistical angle 6 is small, anyons can be described as bosons with a 
hard-core repulsion. Such a system is known to become a superfiuid at low temper- 
atures P[- Similarly, anyons with ^ ~ 27r behave as fermions with a weak, attractive 
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potential. This gives rise to a pairing force between the fermions so they effec- 
tively become bosons which condense. If the particles are charged, we then have a 
superconducting system. These heuristic arguments are actually supported by ap- 
proximative and non-perturbative calculations of anyons at very low temperatures 
0, H and are probably the reason for the breakdown of the perturbative equation 
of state in the same region. 



6 Discussion and conclusion 



Until the discovery of fractional statistics, the thermodynamics of ideal quantum 
gases in two dimensions was completely understood although a few detailed prop- 
erties have first been unravelled during the last couple of years. But these gases are 
in general characterized by a statistical angle which describe anyons interpolating 
between ordinary bosons and fermions. With the recent calculations of the first few 
virial coefficients we have today a fairly good knowledge of the high-temperature 
behaviour of the ideal anyon gas. There are also many indications and arguments 
for the gas undergoing a transition into a superfluid phase at low temperatures. 
However, the details of this transition and the physical properties of the new phase 
are still to a large extent unknown. 

Anyons have many very beautiful properties and there is every reason to be- 
lieve that in the coming years many of the present problems will be much better 
understood or even solved. These results for two-dimensional quantum gases might 
even some time in the future have practical consequences with the ever increasing 
proliferation of new materials and structures in solid state physics. 

We would like to to thank Professor S. Ouvry for several useful suggestions and 
comments. 
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7 Appendix 



For fermions at zero temperature the chemical potential is just the Fermi energy ep- 
As in three dimensions it follows directly form the particle density 



m 



P 



27ih^ Jo 



de 



which gives Sp = 2Ttph^ /m. The corresponding energy density is then 



m 1 



Tip 

m 



(39) 



which also will be the pressure in the fermion gas at zero temperature. 



At temperatures just above zero we can use the Sommerfeld expansion In 
the calculation of the energy, we then find it to terminate after just two terms. 



m 
2^ 



TT 

"6 



(40) 



Similarly, we find no power corrections to the chemical potential, i.e. p = ep as long 
as ksT <^ Ep. This is also in agreement with the exact result (0) which gives 



p = ep + kBT In {l-e-P^'^ 
For the low-temperature energy density of the fermion gas we thus have 

.2/1 \21 



Sp — So 



1 



TT 

y 



pA2 



(41) 



(42) 



where A = (27r/5^^/m)^''^ is the thermal wavelength. Since we have from ([T7|) that 
£p — £b = So, we see that the last term in (^) is the bosonic energy density Sb at 
low temperatures. 

There are no corrections to these results involving higher orders in 1/ pA^ which 
is small at low temperatures. Instead, there will be exponential corrections of the 
order exp(— pA^). One way to see this, is to use the following property of the 
dilogarithmic function 



71 



Li2(x) + Li2(l — x) = — — lnxln(l — x) 



(43) 
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For the energy density (0) of the boson gas we then have 
1 



= 7^Li2(e-^^^) (44) 



^2 



^0 fn' , 2 



pA^ In (l - e-"^') - Li2 (l - g-''^')] . (45) 



(pA2)2 \^ 3 

This provides an exact and exphcit result for the energy density of a two-dimensional 
gas of free bosons valid at all temperatures. Adding the constant (|39|) , we then get 
the corresponding exact result for the fermion gas. 

We now want to show that for a gas in thermodynamic equilibrium the derivative 
of the density p with respect to the to the fugacity z = exp (Pp) is always positive. 
From the grand canonical partition function 

oo 

E=J2 ^""Zm (46) 

where Zj^ is the canonical partition function for exactly particles, we obtain the 
average number of particles 

^1^'^ 1 oo 

{N) = ''-^ = IT.Nz-Z,. (47) 

The density is p = {N)/V where V is the volume of the system. Taking the 
derivative with repect to the fugacity at constant volume, yields 

az ^ j^^Q ^ 
Combining the two terms, we get 

| = -L^,-^-'z,Z.,(iV2-ArAr'). 
Because of the symmetry between A^ and A^' we can rewrite this as 

^ = -^J:^''^'''ZnZ^'C2[n' + n'']-nn') 

N,N' 

N + N' nr ry / 7, r ,r/\2 



5:^^+^'Z^Z^,(Ar-Ar') 



which is always a positive quantity. 
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To lowest order in perturbation theory we had the result ( ^8|) for the pressure 
in the anyon gas as function of the fugacity. The corresponding density is 



1 1 n 



1-'- 



2z 



TV 1 — Z 



(4J 



The last term is seen to give a violation of the above bound at low temperatures 
where the fugacity approaches one. 
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Figure captions: 

Figure 1: Equations of state for free bosons and fermions in two dimensions as 
functions of oc T/p. 

Figure 2: The exchange of two anyons generates a phase-factor which depends 
on the direction of the corresponding rotation in their center of mass. 

Figure 3: Definition of the angle 0^ in the positions of two anyons. 

Figure 4: Spectrum of two anyons in a harmonic oscillator potential. 

Figure 5: The second a) and third b) virial coefficient as a function of the 
statistical parameter a — Qji:. 

Figure 6: The equation of state to first order in 6* for 6* = 0, 0.05, 0.1, with Q 
increasing from below. The temperature at which the approximation breaks down, 
is seen to increase with Q. 
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